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Abstract 

This Quarterly Report analyzes the deployment maneuver of three-axis verti- 
cal constellations with elastic tethers. The deployment strategy devised in 
Quarterly Report #1 has been improved. Dampers have been added to the system. 
Effective algorithms for damping out the fundamental vibrational modes of the 
system have been implemented. This Quarterly Report also shows simulations of a 
complete deployment and a subsequent station keeping phase of a three-mass 
constellation . 
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Figure Captions 


Figure 2.2.1 


Geometry and Lagrangian Coordinates for a Three-Mass 
Constellation 


Figure 2.2.2 

Figure 2.3.1 

Figure 2.6.1 
(a) - (n) 

Figure 2.7.1 
(a) - (1) 


Schematic Model of the Longitudinal Damper 

Schematic Model of the Rotational/Transverse Damper 

Dynamic Response of a Three-Mass Constellation During 
Deployment 


Dynamic Response of a Three-Mass Constellation During 
Station -Keeping 
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1.0 INTRODUCTION 

This is the second Quarterly Report submitted by SAO under contract NAS8- 
36606 , "Analytical Investigation of the Dynamics of Tethered Constellations in 
Earth Orbit (Phase II)," Dr, Enrico Lorenzini, PI, and covers the period from 22 
June 1985 through 21 September 1985. 
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2.0 TECHNICAL ACTIVITY DURING REPORTING PERIOD AND PROGRAM STATUS 

2 . 1 Introductory Remarks 

The two-dimensional equations of motion for three-mass constellations were 
derived in Quarterly Report #1. In that Quarterly Report, however, the two 
tethers connecting the three masses were assumed to be unstretchable . The con- 
trol law used for deployment was a rate control law without angular feedbacks. 
The deployment maneuver was fast but the control law was unable to damp out 
completely the attitude and transverse oscillations of the constellation. This 
time, on the contrary, the two tethers are assumed to be elastic with generic 
characteristics. The right hand terms of the equations of motion developed in 
Quarterly Report #1 are therefore transformed as shown in the following sec- 
tions . 

Two longitudinal dampers, aligned with the two tethers, have been added to 
the system. Angular feedbacks have been implemented in the rate control law in 
order to damp out the attitude and transverse vibrations of the constellations. 
The following sections deal with the incorporation of the above mentioned addi- 
tional features into the mathematical model . The appropriate control parameters 
are then selected and simulations of a deployment and a station-keeping phase 
are performed. 
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2.2 Updated Mathematical Model With Elastic Tethers And Longitudinal Dampers 


The two-dimensional equations of motion (2.2.15) in Quarterly Report #1 
were derived under the general assumption that the three masses were generically 
located in the orbital plane. The only assumption regarding internal forces was 
that a force (tension) Ti was exchanged between masses mi and m 2 and a similar 
force (tension) T 3 was exchanged between masses m 2 and m 3 . In the above-men- 
tioned Quarterly Report the tethers were assumed to be unstretchable and the 
tether tensions were therefore derived by solving the related equations of mo- 
tion for given values of the parameters on the left hand side of the equations. 
By assuming the tethers to be elastic, the tether tension is computed by the 
elastic stretch times the tether stiffnesses as follows: 


Ti = 



tol) 


T3 = ^ (£3 - £03) (2.2.1) 

In equations (2.2.1) , referring to tether #1 and tether #2 as shown in 
Figure 2.2.1, E is the elastic modulus of the tether material, A is the tether 
cross section, £ is the actual tether length and £© is the natural tether 
length. Note that all the quantities referred to tether #2 have the subscript 3 
for consistency of notation. 

Two longitudinal oscillation dampers (one per tether) have been also added 
to the system in order to damp out the elastic longitudinal oscillations of the 
two tethers . In this model the longitudinal dampers are assumed to be two 
classical spring-dashpot systems with a length gain and a velocity gain K^. 


We can therefore write: 
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Ti — Kdi^di + Kdx^di 

T 3 = Kd3^d3 + Kd 3 ^d 3 ( 2 . 2 . 2 ) 

The tethers* tensions, according to equations (2.2.1), and taking into account 
the lengths of the longitudinal dampers, are also given by: 

Ti = ^ (£i - £di - Id) = 

Ts = ^ (^3-^d3-^c3) = Kt3^t3 (2.2.3) 

where EiAx/£i = Ktx and E 3 A 3 /£a = Kta the tethers* stiffnesses and £tii 

are the elastic stretches of the respective tethers. 

By combining equations (2.2.2) and (2.2.3) we finally get the expressions for 
idi and id 3 as follows: 


^dl = (^1 “ ^dl - ^cl) 

^dl 

^d3 = (^3 — ^d3 — ^c3) 


Kdl a 
K-dX 


Kd3 s 

iS 


(2.2.4) 


Equations (2.2.3) and equations (2.2.4) should be added to equations (2.2.15) of 
Quarterly Report #1 in order to get the complete set of equations of motion for 
the system with elastic tethers and longitudinal dampers. The second order 
derivatives £x and £3 must be expressed explicitly in the last two equations of 
motion in order to make them suitable for numerical integration. If we define 7 
= R 3^3 *“ Rx^x so that 7 = Rs £3 - Rx^x (in order to shorten the notation) the 
final equations are as follows 
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Ti = Kti (^1 “ /di ~ ^ci) 


T3 = Kt 3 (^3 “ ^d3 ~ ^c3) 


idi — — ^ {tx — /di — ^ci) — t 

Kdl 


^d3 = ^ (^3 - ^d3 - ^c3) - ^ «d3 

^d3 Kd3 


( 2 . 2 . 5 ) 


B = — l/[mi (£x + 7^1) ™3 (^3 ” 7^3) + R2 (®i + m3)c2j. 

|mi [2£i (6 - n) (£1 + 7) + 302 £icos d sin 0 (£1 + 27) - 3n2£iR2€ cos (2^)] + 
ni3[2£3 (5 — n) (£3 + 7) + 3Q2£3 cos 9 sin ^ (£3 + 27) + 3 Q 2 £3R2€ cos (2^) ] + 
R3 (roi + 1113) [ 2 ee {9 —0) — 30 ^ e 2 cos 9 sin 9 ]| 


€ = [e 2 + 3Q2 csin2^ ] + mtot/ ‘ (207-3n2/Ycos^sin5 ) - 

e (Ti/£i+T 3/£3) /[R2 (mi+m3) ] 


£1 = (biaaa “ ^23x2)/ (axx^aa ~ aia^ai) 
£2 = (b2axx “ bxa2x) / (aix^aa ~ aia^ai) 

where : 

axi = mx (1-Rx) 

ax2 ~ mxR3 

^21 “ m3Rx 

a22 ~ ni3 (I-R3) 

bx = mx [02 ( 3 cos 2 ^-l) (£x+ 27 ) + (^-D) 2 (£x+ 7 ) + 

3n2R2€cosSsin5 + 2nR2e] + Tx 

bi = m3[n2 ( 3 cos 2 ^-l) (£3-27) + (^-n) 2 (£3+7) - 

3n2R2ecos^sin^-2nR2€] + Ts 


(2,2.6) 
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Equations (2.2.5) together with (2.2.6), and with a suitable expression for the 
tether commanded lengths £ci and fc 3 > are the complete mathematical model of the 
system under investigation. Suitable control algorithms will be developed in 
the next subsections . 

2.3 General Considerations Regarding The Damping Of Attitude Librations And 
Transverse Oscillations 

The overall attitude motion of a three-mass system with the middle mass at 
the system c.m. can be well approximated by the two-mass system's equations of 
motion. For the purpose of deriving a tether control algorithm for damping out 
the attitude librations, the above mentioned approximation is more than satis- 
factory. The simplified and linearized equation of motion for in-plane attitude 
oscillations is as follows: 

- 2{6-a)ii + = o (2. 3.1) 

In equation (2.3.1) the second term is the dissipative one. Therefore the 
energy dissipated per each libration cycle is given by: 

T 

- Cl) e dt (2.3.2) 

where r is the period of the in-plane libration. Our goal is to implement a 
tether control law that makes E^ » 0 . 
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I f we consider that ^ is cosinusoidal (with the assumption of t = 0 when ^ = 0) 
and that Q is constant, a good tether control law appears to be: 


I = to (1 - Ked) 


(2.3.3) 


where the gain must be greater than zero and to is the tether length for 6 = 
0. It is important to note that the implementation of the control law (2.3.3) 
assumes that the value of the angle 6 is provided to the reeling systems of the 
constellation. If we now substitute (2.3.3) into (2.3.2), omitting the second 
order terms , we get : 


Ed 


2£oK0 



2Qdt 



(2.3.4) 


In equation (2.3.4) the dominating term is the first integral under parenthesis 
while the second integral is negligible. The orbital rate, therefore, strongly 
influences the damping of the in*“plane librations. It is interesting to note 
that the trajectory followed by the end masses in a damping cycle is not an 
eight-shaped yo-yo cycle but an S-shaped cycle. In other words the tether is 
shortened during the retrograde part of the libration and is lengthened during 
the prograde part as shown in Figure 2.3.1. The control law expressed by equa- 
tion (2.3.3) can also effectively damp the transverse vibrations of the middle 
mass (coordinate e in Figure 2.3.1) by modifying the previous control law for 
the two tethers as follows: 


= toi (1 - Ke<? - K, e//i) 

= to3 (1 - Ke6 + Ke £/^^) 


(2.3.5) 
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Figure 2.3.1 
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The tether control law, therefore, is sensitive both to the overall attitude of 
the constellation (^) and to the transverse displacement of the middle mass (€) . 
A simplified version of (2.3.5) can be obtained by assuming Kc = so that, 
with reference to Figure 2.3.1, we get: 

= lol + €/fi)] = lox (1 - Ml) 

^3 = [1 - e/ia)] - lo3 (1 - Ma) (2.3.6) 

This simplified version, where knowledge of the values of angles and 6^ is 
required, is the one adopted in the following simulations. 


2.4 Selection Of The Control Parameters 

In subsection 2.2 we derived the analytical expression of a control law 
potentially effective in damping out the in-plane oscillations and the trans- 
verse oscillations of the system. In subsection 2.3 we also dealt with the 
addition of passive longitudinal dampers to the system. In this subsection we 
shall preliminarily compute the values of the gains (control parameters) in the 
damping active control laws and in the passive longitudinal dampers. 

In order to maximize the energy transfer between the tethers* longitudinal 
vibrations and the respective longitudinal passive dampers* oscillations the 
natural frequency of each damper should be equal to the natural frequency of the 
associated tether. Unfortunately the tethers* frequencies change during deploy- 
ment because the tethers* lengths change. Since a passive damper is not capable 
of variable tuning, it must be designed for specific tuning. By actively damp- 
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ing longitudinal oscillations by means of reel -control this problem is overcome 
at the expense of further sophistication of the system. This possibility has 
not been eliminated and will be explored in the follow-up to this study. 

Each passive damper has been tuned to the frequency of the associated 
tether at the tether length pertaining to the station -keeping phase. In this 
way the most effective damping of longitudinal vibrations is obtained during 
station-keeping while the damping performance is degraded during deployment. 
The stiffnesses of tether #1 and tether #2 for two millimeter diameter kevlar 
tethers and fully deployed tether lengths are: 

Kti = EiAi/£isk = 67.81 N/m 

Kt 3 = EaAa/^ask = 6.781 N/m (2.4.1) 

In order to tune the longitudinal dampers we must select the dampers* stif- 
fnesses in equations (2.2.2) as follows: 

Kdi = Kti 

Kaa = Kt3 (2.4.2) 


The spring-mass mode longitudinal oscillation frequencies for the constellation 
under investigation in the fully deployed configuration are: 

ojix = \/EiAi/ (£isknii) = 2.735 x 10"^ ra d/ sec 


w^3 = 


x/E^A^/TTa^kTOa) = 2.735 x 10'^ ra d/ sec 


(2.4.3) 
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Some test simulations have been performed in order to identify good damping 
coefficients in ecjuations (2.2.2). By increasing the damping coefficients the 
dampers proved to be more effective. An adimensional coefficient f = 0.9 was 
adopted for both longitudinal dampers although the dynamic response for super- 
critical damping coefficients has not been explored yet. Consequently the di- 
mensional damping coefficients are as follows: 

Kdi = = 4460.24 N/m/ sec 

Kd 3 = 2 ^(jJi 3 m 3 = 446.024 N/m/ sec (2.4.4) 

The selection of the gains (equations 2.3.6) for the rotational/transverse os- 
cillations dampers (active damping system) has been made on a trial and error 
basis. The major constraint is the tether length variation over a damping 
cycle. The maximum tether length variation was finally selected as 1% of the 
fully deployed tether length per degree of system’s libration. The gain in 
equations (2.3.6) is therefore: 


Ke = 0.55 (2.4.5) 

By adopting the above mentioned value the damping of the constellation’s libra- 
tions is very good. A more than satisfactory damping of the system’s transverse 
oscillations is also achieved . The transverse oscillations ’ damping can be 
further improved by using equations (2.3.5) instead of (2.3.6), namely by adopt- 
ing different values for the gains K$ and . This option, however, has not 
been explored yet. Results of the simulation runs performed will be illustrated 


in subsections 2.6 and 2.7. 
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2.5 New Deployment Control Law 

The tether rate control law adopted in the deployment simulations of Quar- 
terly Report #1 has been modified as follows. The initial part of the deploy- 
ment is still an exponentially accelerated phase followed by an exponentially 
decelerated phase with continuous tether speed at the transition between the two 
phases. In formulae we can write: 

Phase I (acceleration) 

^ic = ^11 if fii < fic < £it (2.5.1) 

£ac = £31 simultaneous transition for the two tethers 

Phase II (deceleration) 

£ic ~ (£it “ £if)e"^^ + £if if £ic > £it 

£30 = (£3T - £3f)e-^ + £ 3 , (2.5.2) 

where p = ot/(£if/£iT “ !•) and a = ^ 0sin(2^c) • 

All the characteristic tether lengths pertaining to tether #1 and tether #2, in 
formulae (2.5.1) and (2.5.2), are in the ratio £i.k/£ 3 *k* The angle Oc, is the 
constant in -plane angle assumed by the constellation during the acceleration 
phase. In the deployment control law the deceleration phase does not make the 
tether velocity approach zero. The tethers* controllers simultaneously activate 
the rotational damper when the tether speed in each tether, during the deceler- 
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ating phase, is equal (or very close) to the tether speed required by the rota- 
tional damper for that tether length and that libration angle . The equality of 
the tether speeds, however, does not guarantee that the tether lengths, required 
by the two different control laws, match. For this reason provisions must be 
made for providing a smooth transition . Appropriate values of 6 ^ and of the 
transition length £ix must be adopted in order to reduce the tether length 
mismatch at the transition. In addition a transition control law must be de- 
vised in order to compensate for the residual tether length mismatch. The 
transition control law must be devised in such a way that its initial and final 
velocity are zero. Therefore a semi-cycle of a cosinusoidal law has been adopt- 
ed. All the above mentioned considerations can be translated into formulae as 
follows (similar formulae apply to tether #2) : 

If: [^Xclacc.phmse — [^Iclrot.dmxtp.on that implies -)9(£ic-£if) 


then : 


£lc = £l8k (1 “• ftr ~ KfiS) 


(2.5.3) 


ftr 




(2.5.4) 


and: 


Fotr — (^Ic rot daap.on £lc acc.phasa) /^Isk 


In equation (2.5.4) AT is the transition time. 
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The control parameters which proved to provide a good dynamic response for 
the deployment maneuver are as follows: 

^iT “ 500 m 

= 30° (2.5.5) 

AT = 2000 sec 

A detailed description of the deployment maneuver obtained by using the above 
mentioned control law and control parameters is given in the next subsection. 


2.6 Computer Simulation Of The Deployment Phase 

Several simulation runs of the deployment maneuver have been performed in 
order to select the appropriate control parameters previously described in sub- 
sections 2.4 and 2.5. In this subsection the dynamic response of the system 
during deployment is shown for that "best" selection of parameters. The deploy- 
ment maneuver was started from initial tethers* lengths of in = 20 m and £31 = 
200 m. The system with elastic tether shows a strong tendency to go slack in 
the very beginning of the deployment maneuver . Any small initial mismatch may 
cause a temporary slack tether condition. An in-line thruster and an appropri- 
ately tuned longitudinal damper can help considerably in relieving the slack 
tether tendency. A detailed analysis of the initial deployment phase should 
follow the present study of the deployment maneuver. This type of analysis has 
been postponed, however, for follow-up studies of the three-mass constellations. 
Figure 2.6.1 from (a) to (n) shows the dynamic response during deployment with 
longitudinal dampers tuned to the station-keeping tether lengths and rota- 
tional/transverse dampers activated as described in the previous subsection. 
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Figure 2.6.1 (e) and (f) 
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Figure 2.6.1 (m) and (n) 
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The initial conditions of this deployment maneuver are an in-plane attitude of 
30® with respect to the local vertical and an initial alignment error of the 
three masses (c) of 0,05 m . The dynamic response , however , does not differ 
significantly from a qualitative standpoint for varying values of the above 
mentioned parameters. Figure (a) and Figure (b) show the tether length and 
tether speed respectively of tether #1 vs. time. The same quantities for tether 
#2 can be obtained from those of tether #1 by scaling them up by a factor of 10. 
The deployment is performed in approximately 3 hours. This value however is 
affected by the initial tether length and is therefore eventually affected by 
the position of the reeling system on the Space Station. More important is the 
initial tether speed which must be as close as possible to the initial design 
speed, according to the law £i = a £j , in order to avoid tether slackening. In 
Figure (b) the different phases of the deployment control law are evident: the 
rotational damper-on results in the ripple on the plot at approximately 5500 
sec, while the transition law-off results in the second ripple at 7500 sec. 
Figure (c) and Figure (d) show the tether #1 elastic stretch and longitudinal 
damper's length respectively vs. time. Here again it must be pointed out that 
the system should start the deployment with its natural tether stretch: an 
overstretch or an understretch induces longitudinal oscillations and reduces the 
margin for positive tension in the tether . It is also interesting to note the 
strong coupling between the elastic stretch and the longitudinal damper's 
length, as well as the capability of the system to damp out longitudinal oscil- 
1 at ions as shown by the ripple around 8000 sec . This capabi 1 ity wi 1 1 be shown 
more clearly in the next subsection. Figure (e) shows the in-plane angle (^) 
vs. time. Note that the acceleration phase of deployment is designed for a 
constant in-plane angle. Figure (f) is the phase plane and clearly shows 
the effectiveness of the rotational damper in damping out the attitude libration 
of the constellation. Figure (g) shows the tether tension for tether #1 which 
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is approximately equal to that for tether #2. The tension is proportional to 
the tether elastic stretch (Figure c) . Figure (h) shows the side-view of the 
trajectory of mass #1 (Space Station) with respect to the system center of mass. 
Mass #3 (the ballast) follows a specular trajectory scaled up by a factor 10. 
Mass #2 (the micro-g facility) remains very close to the system c.m. throughout 
the deployment maneuver. Figure (i) is the lateral displacement (e) of mass #2 
vs. time. When the rotational/transverse damper goes on at 5500 sec this oscil- 
lation begins to be damped out. Figure (£) is the phase plane c-€. The damping 
of € is less effective than that of the in-plane angle 6 because the rota- 
tional/transverse damper is tuned to the in-plane oscillation. By using a con- 
trol law like (2.3.5) instead of the simplified (2.3.6), that is, by using a 
multi- frequency damping technique, the damping of the transverse oscillation can 
be further improved. Finally Figure (m) and (n) give the horizontal (along the 
local horizon) and vertical components respectively of the acceleration at the 
micro-g facility vs. time. It is clear from these figures that the initial 
values of the acceleration due to the transverse oscillations of the system and 
the transient oscillations due to snapping of the tethers during transitions in 
the deployment control law are damped out very effectively by the various damp- 
ers. We shall return to this topic again in the next subsection which deals 
with modal vibrations damping during station-keeping. 
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2.7 Computer Simulation Of The Station-Keeping Phase 

The effectiveness of the dampers during station-keeping is shown in the 
following set of plots. These plots have been obtained by simulating the dy- 
namic response of the constellation during a station-keeping phase under the 
following initial conditions. The simulation begins with a tether length for 
tether #1 equal to 909 m and for tether #2 equal to 9090 m; the initial in- 
plane angle is 1® while the initial lateral deflection of mass #2 is 0.10 m. 
Figure 2.7.1(a) shows the tether length for tether #1 vs, time while Figure (b) 
is the side view of mass #1 trajectory. Figure (c) shows the tether’s elastic 
stretch vs. time. Both these quantities are referred to tether #1: tether #2 
has a qualitatively similar behavior. Figure (e) shows the in-plane angle vs. 
time and Figure (f) is the phase plane 6 - 6 , Figure (g) shows the lateral 
deflection (e) of mass #2 (the micro-g platform) vs. time and Figure (h) is the 
phase plane €-€. The saw-teeth- like shape of this plot, near the beginning of 
the simulation, is due to an unsufficiently small plotting step and therefore it 
has no dynamic meaning. It is evident from these plots that any initial pertur- 
bations of the constellation are abated very effectively by the dampers. The 
final result is shown in Figure (i) and (1) which give the horizontal and verti- 
cal component respectively of the acceleration at mass #2 vs. time. The maximum 
accelerations at the end of this simulation are much smaller than 10'®g. This 
last point must not be misinterpreted: it means only that the system can effec- 
tively abate transient vibrations while the steady state acceleration level at 
mass #2 will depend upon the external steady state perturbations. These pertur- 
bations are not simulated in this analysis which deals with the transient re- 
sponse of the three-mass constellation. 
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Figure 2.7.1 (c) and (d) 
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2.8 Concluding Remarks 

Passive dampers for tethers* longitudinal oscillations have been added to 
the system with elastic tethers. Active control of rotational/transverse oscil- 
lations have been also added to the constellation . A new deployment strategy 
which drives the system to a steady state condition in approximately 3 hours has 
been devised. The initial part of deployment with elastic tethers, however, 
requires a more detailed investigation; as a result of such investigation the 
deployment duration may increase depending upon the minimum separation length of 
the end masses at the start of deployment. The dampers have been designed to 
provide an effective damping of the longitudinal, libratory and transverse os- 
cillations of the system. The effectiveness of the dampers has been proven by 
computer simulations of the deployment and station-keeping phase. The damping 
of the modal vibrations can be further improved by multi- frequency active damp- 
ing. 


3.0 PROBLEMS ENCOUNTERED DURING REPORTING PERIOD 


None 
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4.0 ACTIVITY PLANNED FOR THE NEXT REPORTING PERIOD 


During the next reporting period we will analyze the dynamics of three-mass 
constellations when the middle mass travels along the tether . Stability con- 
straints and appropriate control laws will be sought. The investigation of the 
dynamics of a five-mass system, vertically oriented, attached to the Shuttle 
will also be initiated. This system potentially could be used as a scientific 
platform for measuring geophysical gradients. 



